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Abstract 

Using two alternative [0/H]-[Fe/H] dependences, the empirical oxygen 
abundance distribution (EGD) is deduced from two different samples 
involving (i) 268 K-giant bulge stars (Sadler et al. 1996), and (ii) 
149 globular clusters (Mackey & van den Bergh 2005) for which the 
iron abundance distribution is known, in addition to previous results 
(Caimmi 2001) related to (hi) 372 solar neighbourhood halo subdwarfs 
(Ryan & Norris 1991). Under the assumption that each distribution 
is typical for the corresponding subsystem, the EGD of the Galactic 
spheroid is determined weighting by mass. The trend is non mono- 
tonic, with the occurrence of one minimum and two maxima within 
the domain. The data are fitted, to an acceptable extent, by simple 
models of chemical evolution implying both homogeneous and inho- 
mogeneous mixing, provided star formation is inhibited during halo 
formation and enhanced during bulge formation, with respect to the 
disk. The theoretical oxygen abundance distribution (TGD) is first 
determined for the halo and the bulge separately, and then for the 
Galactic spheroid weighting by mass. Though a G-dwarf problem 
seems to exist for both the halo and the bulge, it could be allevi- 
ated or removed allowing an increasing star formation efficiency dur- 
ing the formation of a Galactic subsystem, which does not necessarily 
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imply gas infall. The results are independent of the power-law ini- 
tial mass function (IMF) exponent, provided both the lower stellar 
mass limit and the mass fraction in long-lived stars and remnants per 
star generation, are suitably changed. Then the effect of star forma- 
tion inhibiting or enhancing gas is counter-balanced. Simple models 
implying homogeneous mixing are unable in fitting the empirical age- 
metallicity relation (EAMR) recently determined from a homogeneous 
sample of globular clusters (De Angeli et al. 2005), which shows a 
non monotonic trend characterized by large dispersion. On the other 
hand, simple models involving inhomogeneous mixing yield a theoret- 
ical age-metallicity relation (TAMR) which reproduces the data to an 
acceptable extent. With regard to gas ouflow from the proto-halo, 
acceptable models make different predictions according if the Galactic 
spheroid and disk underwent separate or common chemical evolution. 
In the former alternative, less than one third of the bulge mass out- 
flowed from the proto-halo. In the latter alternative, the existence 
of an unseen baryonic halo (or equivalent amount of gas lost by the 
Galaxy) with mass comparable to bulge mass, is necessarily needed. 
In this view, the outflowing proto-halo gas which remains bound to 
the Galaxy, makes both the bulge and the disk. The treshold star 
mass below which the halo is not visible (or an equivalent amount of 
gas has been lost) is calculated as tjiq rs 0.25 m for IMF exponent 
p = 2.9 and mo ~ 0.10 m Q for p = 2.35. Conversely, p 2.8 for lower 
limit stellar mass, m m f = 0.25 m and p « 2.6 for niQ ss 0.10 m Q . 
keywords - galaxies: evolution - stars: formation; evolution. 

1 Introduction 

The existence of a G-dwarf problem i.e. the observation of too few metal de- 
ficient G dwarfs (or, more generally, of a selected spectral type) with respect 
to what expected from the Simple model of chemical evolution (e.g., Searle 
& Sargent 1972; Pagel & Patchett 1975; Haywood 2001) was first established 
in the solar neighbourhood (van den Bergh 1962; Schmidt 1963). Though 
in less extreme form, a G-dwarf problem appears to exist in both the halo 
(e.g., Hartwick 1976; Prantzos 2003) and in the bulge (e.g., Ferreras, Wyse & 
Silk 2003). In addition, a G-dwarf problem has been detected in both bulge- 
dominated and disk-dominated galaxies (Henry & Worthey 1999), which is 
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consistent with the idea that the G-dwarf problem is universal (Worthey, 
Dorman & Jones 1996). 

The deficit of metal-poor stars (with respect to the prediction of the 
Simple model) may be interpreted in different ways, such as changes in the 
initial mass function (Schmidt 1963; Adams & Fatuzzo 1996; Bromm 2004; 
Bromm & Larson 2004; Larson 2005), inflow of unprocessed (Larson 1974) 
or processed (Thacker, Scannapieco & Davis 2002) material from outside, 
or evolution with inhomogeneous mixing (Searle 1972; Malinie et al. 1993; 
Caimmi 2000, 2001b, hereafter quoted as C00 1 and C01, respectively; Oey 
2003; Karlsson 2005). For additional alternatives and further details see e.g., 
Pagel & Patchett (1975); Pagel (1989). 

In addition to the G-dwarf problem, a lack of a well-defined empirical 
age-metallicity relation (EAMR) seems to be established for both the disk 
solar neighbourhood (e.g., Meusinger, Reimann & Stecklum 1991; Edvards- 
son et al. 1993; Rocha-Pinto et al. 2000; Feltzing, Holmberg & Hurley 2001; 
Nordstrom et al. 2004; Karatas, Bilitz & Schuster 2005) and globular clusters 
(Salaris & Weiss 2002; De Angeli et al. 2005). The large scatter observed 
in the EAMR is probably universal, at least with regard to massive enough 
(M ~ lO lo m ) galaxies, independent of the morphological type. 

Inhomogeneous (i.e. implying inhomogeneous gas mixing) models of 
chemical evolution succeed in both providing a solution to the G-dwarf prob- 
lem and reproducing substantial scatter exhibited by the EAMR. The current 
paper aims to investigate if inhomogeneous simple models of chemical evo- 
lution are also consistent with the metallicity distribution in the Galactic 
spheroid, deduced weighting by mass data belonging to subsystems, more 
specifically solar neighbourhood halo subdwarfs, K-giant bulge stars, and 
globular clusters. In this view, the next step is to see what constraints are 
related to the formation and the evolution of the Galaxy. 

With regard to halo subdwarfs, the oxygen abundance has been deduced 
(C01) from data related to a sample of 372 kinematically selected halo stars 
(Ryan & Norris 1991) after conversion of [Fe/H] into [O/H] using two alterna- 
tive empirical relations, involving the presence or absence of [O/Fe] plateau 
for sufficiently low [Fe/H] values. 

1 With regard to this reference, two points may carefully be kept in mind, namely 
(i) values of a few parameters must be corrected as explained in Caimmi (2001b), Sect. 3, 
second paragraph therein, and (ii) the majority of figures do not correspond to the caption, 
as explained in the erratum (Caimmi 2001a). 
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With regard to globular clusters and bulge stars, a similar procedure is 
applied to: (i) a homogeneous sample of 55 objects for which the EAMR is 
also known (De Angeli et al. 2005), (ii) an inhomogeneous sample of 149 
objects (Mackey & van den Bergh 2005), and (iii) a homogeneous sample of 
268 bulge K-giant (Sadler, Rich & Terndrup 1996), as reported in Section 2. 
In addition, it is derived therein a putative oxygen abundance distribution in 
all halo objects, under the assumption that both globular clusters and field 
stars underwent a common chemical evolution. The same is made also for 
the Galactic spheroid, with the inclusion of bulge stars. A comparison with 
the predictions of simple models, involving homogeneous and inhomogeneous 
mixing, is made in Section 3 and 4, respectively. The discussion and the 
conclusion are the subject of Section 5 and 6, respectively. 

2 The data 

2.1 Empirical age-metallicity relation in globular clus- 
ters 

Accurate relative ages for a sample of 55 Galactic globular clusters have 
recently been determined (De Angeli et al. 2005). The ages were obtained 
by measuring the difference between the horizontal branch and the turnoff 
in two internally photometrically homogeneous databases with 16 objects in 
common. The related EAMR is derived from absolute ages and related errors 
(De Angeli 2005) for objects belonging to a single database. With regard to 
objects in common where a value and an error exist for each database, the 
absolute age is calculated as the centre of the intersection of the two intervals, 
and the error as the corresponding semiamplitude. For instance, 11.42=1=0.34 
and 10.96 =F 0.60 yield 11.32 =f 0.24 2 . Following De Angeli et al. (2005), 
metallicities are calibrated over two different scales, namely CG (Carretta & 
Gratton 1997, as extended by Carretta et al. 2001) and ZW (Zinn & West 
1984). In addition to the above homogeneous values, a third alternative 
consists in using the data (taken from different sources) from Harris (1996) 
catalogue (2003 update), also listed (with one exception) in a recent paper 
(Mackey & van den Bergh 2005), where five classification types are defined. 

2 The absolute age could be calculated as a weighted mean, provided the errors available 
for each database may be related to empirical variances. 
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With regard to the current sample of 55 objects, three classification types 
among the above mentioned five shall be retained, namely: BD for bulge/ disk 
clusters (5 objects); OH for old halo clusters (36 objects); and YH for young 
halo clusters (13 objects). For further details, see Mackey & van den Bergh 
(2005). An additional cluster (NGC 6366) has to be considered by itself, as 
belonging to BD type ([Fe/H]> —0.8) with respect to CG and ZW metallicity 
scales, and to OH type ([Fe/H]< —0.8) with resepct to the value listed in 
Harris catalogue (Mackey & van den Bergh 2005). 

The EAMR for the sample under discussion is represented in Fig. 1, with 
regard to the above mentioned metallicity scales. An inspection to Fig. 1 
discloses the following. 

(i) Age differences related to different metallicity scales are small. 

(ii) The formation of globular clusters in the Galactic halo was a continuous 

process which span over 4 Gyr. Low-metallicity halo clusters were 
generated over m 2 Gyr. Intermediate-metallicity halo clusters were 
generated over 4 Gyr. High-metallicity halo clusters and bulge/disk 
clusters were generated over ps 1 Gyr. 

(iii) No clear distinction between "old halo" and "young halo" clusters seems 

to be in terms of age difference. Perhaps it should be better referring 
to "primeval halo" and "accreted halo" clusters. The inclusion of YH 
clusters has no appreciable effect on the EAMR. 

Accordingly, the following picture of Galactic evolution may be inferred. 
Low-metallicity stars formed within the proto-halo and after about 1 Gyr 
within the proto-bulge/proto-disk. Bulge formation lasted about 1 Gyr, con- 
sistent with the absence of bulge stars younger than 10 Gyr (Zoccali et al. 
2003), while halo formation ended about 8 Gyr ago, similar to what has been 
found for the thick disk (Liu & Chaboyer 2000; Ibukiyaua & Arimoto 2002). 
Then halo and thick disk formation lasted about 4 Gyr and 3 Gyr, respec- 
tively. In addition, one has to keep in mind that different specific angular 
momentum distributions occur in halo-bulge and thick disk-thin disk stars, 
which implies that the two subsystems had different dynamical evolutions 
(Ibata & Gilmore 1995). 
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2.2 Metallicity distribution in globular clusters and 
field stars 



The empirical distribution of oxygen abundance in globular clusters is de- 
duced from a sample of 149 globular clusters where [Fe/H] has been deter- 
mined from different sources (Mackey & van den Bergh 2005) and from a 
sample of 55 globular clusters where [Fe/H] has been determined from a sin- 
gle source (De Angeli et al. 2005). In dealing with simple models of chemical 
evolution, involving the assumption of instantaneous recycling, the predicted 
metal abundance has to be compared with the observed oxygen abundance 
(e.g., Pagel 1989; COO; C01). Unfortunately, oxygen is more difficult than 
iron to detect, and an empirical relation is needed, to express the former as 
a function of the latter. 

To this respect, a clear dichotomy appears to be among authors who sup- 
port a plateau in [O/Fe] for stars with [Fe/H]~ —1 (e.g., Carretta, Gratton & 
Sneden 2000); and those who do not (e.g., Israelian et al. 2001); to get further 
insight, see the proceedings edited by Barbuy et al. (2001). For exploiting 
both the above possibilities into consideration, the following relations (C01) 
shall be used: 
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in presence or in absence of [Fe/H] plateau, respectively, where the oxygen 
solar abundance is taken to be O = 0.0056 (Allende-Prieto, Lambert, & 
Asplund 2001); for further details, see C01. 

A number of [0/H]-[Fe/H] relations lying between those expressed by 
Eqs. (1) and (2), have been derived from recent investigations (Jonsell et 
al. 2005; Fulbright, Rich & McWilliam 2005; Garcia Perez et al. 2005; 
Melendez et al. 2005). To allow comparison with previous results related to 
halo subdwarfs (C01), Eqs. (1) and (2) shall be used in the current attempt, 
and hereafter quoted as "in presence" and "in absence" of [O/Fe] plateau, 
respectively, with regard to sufficiently low metallicities, [Fe/H]~ —1. 

On the other hand, extremely metal deficient ([Fe/H]< —4) stars are 
known to be oxygen overabundant (e.g., Christlieb et al. 2002; Iwamoto et 
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al. 2005), [0/Fe]« 2.5 (Bromm & Loeb 2003; Frabel et al. 2006). Then 
Eqs. (1) and (2) cannot be used in this case, and the chemical evolution shall 
be restricted to a later epoch where only pop. II star formation occurred, 
[Fe/H]> -4, say. 

Let us define the oxygen abundance normalized to the solar value, 0, as: 



log = log 



_0_ 
0~ 







(3) 



and let Alog0 = A[0/H] = [0/H]+- [O/H] - be a logarithmic, oxygen abun- 
dance bin deduced from A[Fe/H] by use of Eq. (1) or (2). The related, oxygen 
abundance bin is: 



A0 = A + + A~0; A T 0=|0-0 T | 



T = exp 10 



H 



(4a) 
(4b) 



where in general, exp^ defines the power of basis £ and, in particular, exp 
defines the power of basis e, according to the standard notation. As in C01, 
bins in [Fe/H] equal to 0.2 dex shall be used (e.g., Norris & Ryan 1991; 
Huchra, Brodie & Kent 1991; Perrett et al. 2002). 

The empirical, differential metallicity distribution (hereafter referred to 
as EGD) in a selected class of objects, is defined as (Pagel 1989; COO; C01): 



iP(<j)T A t 0) =log 



AN 

NA(f) 



(5) 



where A0 is the bin width, AiV is the number of sample objects with oxygen 
abundance belonging to a bin centered in 0, and N is the total number of 
sample objects. The differential distribution is used instead of the cumulative 
distribution, as it is a more sensitive test (Pagel 1989) and allows direct 
comparison between different samples. The uncertainty on AiV has been 
evaluated from Poisson errors (e.g., Ryan & Norris 1991), as A(AiV) = 
(AiV) 1//2 , and the related uncertainty in the EGD is (e.g., C01): 
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where ip~ diverges to — oo in the limit AN — > 1. For further details, see C01. 

The [Fe/H]-[0/H] relation and corresponding mean fractional oxygen 
abundance, 0, and half bin width, A T (j), in presence of [O/Fe] plateau (PP), 
according to Eq. (1), and in absence of [O/Fe] plateau (AP), according to 
Eq. (2), respectively, are shown in Tab. 1 for the metallicity range of interest. 
The EGD derived from the sample studied by De Angeli et al. (2005), using 
Eqs. (3), (4), (5), (6), is listed in Tab. 2 in presence of [O/Fe] plateau with 
regard to CG metallicity calibration, and in absence of [O/Fe] plateau with 
regard to ZW metallicity calibration. 

The related plots are shown in Fig. 2, top left and top right, respec- 
tively. The dotted vertical line marks the transition from halo (OH, YH) to 
bulge/disk (BD) morphological type. Bottom left and bottom right panels 
represent a reduced sample of 42 objects, where YH clusters have been re- 
moved. The inclusion of YH clusters (safely stripped from accreted dwarf 
galaxies like Sagittarius) appears to have no appreciable effect on the EGD. 
The distribution is bimodal with the occurrence of two maxima, close to 
the beginning of the evolution and to the transition from halo to bulge/disk 
morphological type, respectively. 

The EGD derived from the sample studied by Mackey & van den Bergh 
(2005), using Eqs. (3), (4), (5), (6), is listed in Tab. 3 in presence and in ab- 
sence of [O/Fe] plateau, with regard to metallicity values taken from different 
sources. 

The related plots are shown in Fig. 3, top left and top right, respectively. 
Bottom left and bottom right panels represent a reduced sample of 107 ob- 
jects, where only OH clusters have been retained. The inclusion of YH clus- 
ters appears to have no appreciable effect on the EGD which, on the other 
hand, is bimodal with the occurrence of two maxima, close to the beginning 
of the evolution and to the transition from halo to bulge/disk morphological 
type, respectively. 

The EGD derived from the sample of 268 K-giant bulge stars in Baade's 
Window studied by Sadler, Rich & Terndrup (1996) 3 , using Eqs. (3), (4), (5), 
(6), is listed in Tab. 4 both in presence and in absence of [O/Fe] plateau. 

The related plots are shown in Fig. 4, bottom left and bottom right, re- 
spectively. Bottom panels represent the EGD deduced from a sample of 372 

3 The abundance distribution of the sample stars has recently been revised, but with 
no substantial change (Fulbright ct al. 2005, 2006). 
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Table 1: The [Fe/H]-[0/H] relation and corresponding mean fractional oxy- 
gen abundance, 0, and half bin width, A T <fi, in presence (PP) and in absence 
(AP) of [O/Fe] plateau. To save space, F stays for [Fe/H] and O for 3[0/H]. 
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Table 2: The empirical, differential metallicity distribution (EGD) in glob- 
ular clusters, deduced from a sample of 55 objects studied by De Angeli et 
al. (2005), in presence of [O/Fe] plateau (PP) with regard to CG metallic- 
ity calibration, and in absence of [O/Fe] plateau (AP) with regard to ZW 
metallicity calibration. 
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Table 3: The empirical, differential metallicity distribution (EGD) in globular 
clusters, deduced from a sample of 149 objects studied by Mackey & van den 
Bergh (2005), both in presence (PP) and in absence (AP) of [O/Fe] plateau. 
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Table 4: The empirical, differential metallicity distribution (EGD) in K-giant 
bulge stars in the Baade's Window, deduced from a sample of 268 objects 
studied by Sadler et al. (1996), both in presence (PP) and in absence (AP) 
of [O/Fe] plateau. 
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1. 


,604 


-0 


.444 


0.069 


0.060 


46 


1 


.129 


-0 


,275 


1. 


.180 


-0 


.417 


0.079 


0.067 


36 


0. 


.897 


-0 


,175 


0. 


,868 


-0 


.284 


0.079 


0.067 


36 


0. 


.713 


-0 


.169 


0. 


.638 


-0 


.245 


0.089 


0.074 


29 


0. 


.566 


-0 


,100 


0. 


.470 


-0 


.142 


0.093 


0.076 


27 





.450 


-0. 


.201 


0, 


,341 


-0 


.210 


0.121 


0.094 


17 





.357 


-0 


.377 


0, 


,254 


-0 


.353 


0.176 


0.125 


9 





.284 


-0. 


.930 


0, 


,187 


-0 


.873 


0.533 


0.232 


2 





.205 


-1. 


.385 


0. 


.137 


-1 


.040 


oo 


0.301 


1 





.129 


-1. 


.185 


0, 


,101 


-0 


.907 


oo 


0.301 


1 


0. 


.081 





.000 


0. 


.074 





.000 


0.000 


0.000 








.051 


-0. 


.785 


0, 


,055 


-0 


.640 


oo 


0.301 


1 
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halo subdwarfs in the solar neighbourhood, studied by Ryan & Norris (1991); 
for further details, see C01. 



3 Inferred metallicity distribution in the Galac- 
tic spheroid 

The Galactic spheroid may safely be conceived as made of three main sub- 
systems, namely: (i) globular clusters, (ii) field halo stars, and (iii) field 
bulge stars. Let us assume a bulge mass, Mb = lO lo m (e.g., Wyse & 
Gilmore 1992; Kent 1992; Prantzos & Silk 1998), a (baryonic) halo mass, 
Mu = 10 9 m o (e.g., Carney, Latham & Laird 1990; Wyse & Gilmore 1992), 
and a ratio of globular clusters to halo mass, M G c/M H = 1/100 (e.g., Li & 
Burstein 2003), where the above values of bulge and halo mass may safely be 
used including or not globular clusters. Accordingly, the mass of the Galactic 
spheroid and the fractional mass of related subsystems read: 

M = M GC + M FH + M FB = (0.001 + 0.1 + 1) 10 10 m o = 1.101 m Q ; (7a) 

(7b) 

(7c) 

(7d) 

where the indices, GC, FH, FB, mean globular clusters, field halo stars, and 
field bulge stars, respectively. 

Let N be the total number of long-lived (i.e. life time longer than the age 
of the Galaxy) stars in the Galactic spheroid, and AiV the number of long- 
lived stars within a selected metallicity bin. The relative frequency, AN/N, 

AiV _ N GC AN gc iVFH AiVFH N FB AN FB 
N N N GC N N FH N Np-Q ' U 

where ANxy/Nxy is the relative frequency belonging to XY subsystem, XY 
= GC, FH, FB, with regard to the selected metallicity bin. Samples related 
to different subsystems may be (and in fact are) made of different objects, 
provided each sample is statistically significant. 



M GC 


1 


= 0.0009 


M 


1101 


M FH 


100 


= 0.0908 


M 


1101 


M FB 
M 


1000 
1101 


= 0.9083 
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With regard to field (halo and bulge) stars, let us define a selected spectral 
class of long-lived stars as belonging to a selected mass range, m\ <m < m 2 , 
and suppose it does not significantly depend on the evolution, as in disk stars 
(Rocha-Pinto & Maciel 1997). Then the fractional number of stars belonging 
to an assigned mass range, equals the fractional number of stars belonging 
to the whole mass domain, provided the initial mass function (IMF) did not 
significantly change within the Galactic spheroid. Accordingly, the following 
relation holds: 



(%0 = ^ XY = GC, FH, FB ; (9) 

\ -iV /mi ,mo -/V 



'Nxy\ 
and, on the other hand: 

Nxy _ rnxyNxv _ M XY 
N ~ mN ~ M 



(10) 



where m X y an d m are the mean mass of long-lived stars in XY subsystem and 
in the Galactic spheroid, respectively, and Wixy = m owing to the assumption 
of universal IMF. 

The combination of Eqs. (8), (9), and (10) yields: 

AN _ M GG AN GG M FH AiV F H M FB AN FB 
N M N GC M iV F H M A'li! 

and the related uncertainty is obtained using the standard formula of linear 
propagation of errors together with evaluation of Poisson errors (e.g., Ryan 
& Norris 1991), A(AA^xy) = (AiV XY ) 1/2 . The result is: 

A AN _ M GC (AN GC ) 1 / 2 M FH (AA^) 1 / 2 M FB (AN^) 1 / 2 
N M N GC M .V|,„ M N FB ' 

which may explicitly be calculated using Eqs. (7) and the data listed in Tab. 3, 
Tabs. 1-2 in C01, and Tab. 4, concerning globular clusters, field halo stars, 
and field bulge stars, respectively. 

The EGD in stars of the Galactic spheroid results from the combination 
of Eqs. (5), (6), and (11), as: 



ip = log 



M GC AN gc , M FH AiV F H , M FB AAT FB 



M N GC A(j) M .\|„Ao M .\ K |.Ao 



(13a) 
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log 



' MAN/N) 
[ + AN/N j 
T V [ M XY AiV XYT (AiV XY ) 1 / 2 ' 
5 ' M iV XY A0 



(13b) 
(13c) 



where XY = GC, FH, FB. 

As a first application of Eqs. (11), (12), and (13), let us take into consid- 
eration the light Galactic spheroid, made of globular clusters and field halo 
stars. To this aim, the globular cluster sample (Mackey & van den Bergh 
2005) has been reduced to N = 148 with the exclusion of Liller I owing to its 
high metal content, [Fe/H]=0.22, which exceeds the maximum metallicity in 
solar neighbourhood halo subdwarf sample (Ryan & Norris 1991) used here. 
The resulting EGD is listed in Tab. 5 both in presence and in absence of 
[O/Fe] plateau. The related plots are shown in Fig. 5 (top panels), left and 
right, respectively. 

A second application concerns the massive Galactic spheroid, including 
globular clusters, field halo stars, and field bulge stars. The resulting EGD 
is listed in Tab. 6 both in presence and in absence of [O/Fe] plateau. The 
related plots are shown in Fig. 5 (bottom panels), left and right, respectively. 

The EGD in the Galactic spheroid is similar to its counterpart related to 
globular cluster subsystem, Fig. 3, which is bimodal with the occurrence of 
two distinct maxima, and a minimum soon before the transition from halo 
to bulge/disk morphological type, represented by a dotted vertical line. It is 
apparent that (closed or open) simple models of chemical evolution cannot 
provide a satisfactory explanation to a bimodal EGD, and a different model 
is needed. 

Accordingly, let us suppose the halo and the bulge underwent distinct 
chemical evolutions, calculate the related theoretical differential metallicity 
distribution (hereafter referred to as TGD) in a selected spectral class of long- 
lived stars, the resulting TGD, and compare with its empirical counterpart. 
To this aim, simple models implying both homogeneous and inhomogeneous 
mixing shall be used, as in C01. 
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Table 5: The empirical, differential metallicity distribution (EGD) in the 
light Glalctic spheroid (globular clusters and field halo stars), deduced from 
a sample of 148 globular clusters (Mackey & van den Bergh 2005) and a 
sample of 372 field halo subdwarfs (Ryan & Norris 1991), both in presence 
(PP) and in absence (AP) of [O/Fe] plateau. 



PP AP 










i> 











A--0 






\N 
N 


A— 

^ N 





.958 


-2. 


.048 





.951 


-2, 


.167 


1.336 


0.291 





.003 


0.003 


0. 


.713 


-1. 


.080 


0. 


.638 


-1, 


.156 


0.258 


0.161 





.014 


0.006 





.566 


-0 


.821 





.470 


-0 


.864 


0.202 


0.137 





.020 


0.007 





.450 


-0. 


.496 





.345 


-0 


.505 


0.148 


0.110 


0, 


.033 


0.009 





.357 


-0. 


.367 





,254 


-0 


.342 


0.142 


0.107 





.035 


0.010 





.284 





.062 





,187 


0, 


.120 


0.092 


0.076 


0, 


,075 


0.014 





.205 





.096 





,137 


0, 


.441 


0.072 


0.062 


0, 


,116 


0.018 


0. 


.129 





.371 





,101 


0, 


.650 


0.066 


0.057 


0, 


,137 


0.019 





.081 





.653 





,074 





.865 


0.059 


0.052 


0, 


.166 


0.021 





.051 





.707 





,055 





.852 


0.071 


0.061 


0, 


.118 


0.018 





.032 





.792 





,040 


0, 


.870 


0.082 


0.069 





,091 


0.016 





.020 





.840 





.030 





.852 


0.100 


0.081 


0, 


.064 


0.013 





.013 





.980 





.022 





.925 


0.107 


0.086 


0, 


.056 


0.012 





.008 





.972 


0. 


.016 





.850 


0.141 


0.106 


0, 


.035 


0.010 





.005 





.961 


0. 


.012 





.773 


0.189 


0.131 


0, 


.021 


0.008 





.002 





.700 


0. 


.007 





.385 


0.228 


0.149 


0, 


.016 


0.007 
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Table 6: The empirical, differential metallicity distribution (EGD) in the 
massive Galactic spheroid (globular clusters, field halo stars, and field bulge 
stars), deduced from a sample of 149 globular clusters (Mackey & van den 
Bergh 2005), a sample of 372 field halo subdwarfs (Ryan & Norris 1991), 
and a sample of 268 field K-giant bulge stars (Sadler et al. 1996), both in 
presence (PP) and in absence (AP) of [O/Fe] plateau. 



PP I AP 










1> 









4 


A--0 






AN 

N 




AN 


4. 


.496 


-2 


,172 


7, 


,443 


-2, 


.515 


0, 


533 





oo o 

,232 





,007 


0. 


005 


3 


.572 


-2 


,072 


5, 


,476 


-2, 


.381 


0, 


r - o o 

533 


0. 


oo o 

232 





,007 


0, 


005 


2 


.837 


-1 


,370 


4 


.028 


-1 


.646 


0. 


,189 


0, 


131 


0, 


,028 


0. 


010 


2 


.254 


-0 


.969 


2 


.963 


-1 


.211 


0. 


,125 


0, 


,097 





.055 


0. 


,014 


1 


.790 


-0 


.611 


2 


.180 


-0 


.820 


0. 


089 


0, 


,074 





.100 


0. 


019 


1 


.422 


-0 


.310 


1 


.604 


-0 


.486 


0. 


069 


0, 


060 





.159 


0. 


023 


1 


.129 


-0. 


.316 


1 


,180 


-0 


.458 


0, 


080 


0. 


068 


0, 


.125 


0. 


021 





.897 


-0 


.217 





,868 


-0. 


.326 


0, 


,079 


0. 


,067 


0, 


.125 


0. 


,021 


0. 


.713 


-0. 


.205 





,638 


-0 


.281 


0, 


,091 


0. 


.075 


0, 


.102 


0. 


019 


0. 


.566 


-0. 


.134 





,470 


-0. 


.176 


0, 


,095 


0. 


,078 


0, 


,095 


0. 


,019 





.450 


-0 


.221 





,341 


-0 


.230 


0, 


,122 


0. 


,095 


0, 


,062 


0. 


,015 





.357 


-0 


.376 





,254 


-0. 


.352 


0, 


,173 


0. 


.127 





,034 


0. 


Oil 





.284 


-0 


,673 


0, 


,187 


-0, 


,615 


0, 


.259 


0. 


.161 





,014 


0. 


006 


0. 


.205 


-0. 


.819 





,137 


-0, 


,474 


0, 


,195 


0. 


,134 





,014 


0. 


,005 


0. 


129 


-0 


.561 





.101 


-0. 


,282 


0. 


,171 


0, 


,122 


0, 


.016 


0, 


005 


0. 


.081 


-0 


.384 


0, 


.074 


-0. 
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0. 


,059 


0, 


,052 


0, 


.015 


0, 


002 


0. 


.051 


-0 


.211 


0, 


.055 


-0. 


066 


0. 


,192 


0, 


,132 


0, 


.014 


0. 


005 


0. 


032 


-0 


.245 





.040 


-0. 


,167 


0. 


,082 


0, 


,069 


0, 


.008 


0, 


001 


0. 


020 


-0 


.197 





.030 


-0. 


,185 


0. 


,100 


0, 


,081 


0, 


.006 


0, 


001 


0. 


.013 


-0 


.057 


0, 


.022 


-0, 


,112 


0, 


,107 


0, 


,086 


0, 


,005 


0, 


001 


0. 


.008 


-0 


.065 


0. 


,016 


-0, 


,187 


0, 


,141 


0, 


106 





,003 


0, 


001 


0. 


.005 


-0 


.076 


0, 


.012 


-0, 


,265 


0, 


,189 


0, 


131 


0, 


,002 


0, 


001 


0. 


.002 


-0 


,337 





,007 


-0. 


,653 


0, 


.228 


0. 


,149 


0, 


,001 


0, 


001 
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4 Homogeneous simple models 



With regard to homogeneous simple models with star formation inhibiting 
gas (e.g., Hartwick 1976; COO; C01), the TGD is represented as a straight 
line (e.g., Pagel 1989; COO; C01): 

diV 

^ = log ivd^ = a0 + b ; (14) 

and the explicit expression of the coefficients, a and b, reads: 

1 O 



In 10 p" ' 



(15) 



Ho O 



b = \og -acf> ; (16) 

\fi - n f p" J 

where O is the solar oxygen abundance, p" the effective (oxygen) yield, fi 
the (allowing star formation) gas mass fraction, and the indices, o and /, 
denote the beginning and the end of evolution, respectively. 

The oxygen abundance, O, may be related, to a good extent, to the gas 
mass fraction, //, as (e.g., COO; C01): 

O-O = -j3"ln^ ; (17) 

f = ^- K ; (18) 

where p, is the real (oxygen) yield and k is the ratio of (inhibiting star 
formation) gas mass fraction to long-lived star and remnant mass fraction. 
In addition, the following assumptions have been made: (i) instantaneous 
recycling; (ii) universal power-law IMF; and (iii) gas inhibition from forming 
stars at a rate proportional to the star formation rate. For further details 
see e.g., Hartwick (1976); Pagel (1989); COO; C01. 

A deeper analysis shows that, given a simple model with star formation 
inhibiting gas, long-lived lower stellar mass limit, m m f, and inhibition pa- 
rameter, k, it is equivalent to a simple model with same values of independent 
parameters except the above mentioned two, (m m f) 1 < m m f and K-y — < k; 
and vice versa. Accordingly, the related yield, pi, has the same value as the 
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effective yield, p", expressed by Eq. (18). For an exhaustive discussion, see 
C01. 

Simple models with star formation inhibiting gas (Hartwick 1976) can 
be generalized to negative values of the inhibition parameter, k < 0, where 
star formation is enhanced by additional gas with same composition as the 
pre-existing one, at a rate: 

dD dS 

— = ; «<0; (19) 

which corresponds to a total amount: 

D(t) -D = aK[S(t) - S ] = K [s(t) - s ] ; « < ; (20) 

where D, S, and s represent the mass fraction in additional gas, stars formed, 
and long-lived stars and remnants, respectively, and a is the mass fraction of 
a star generation which remains locked up in long-lived stars and remnants. 

The additional gas may be conceived as inflowing from outside, with 
same composition as the pre-existing gas, and entirely turned into stars. 
Accordingly, the mass fraction with respect to the original system, u, + s, 
reads: 

H(t) + s{t) = 1 - D(t) ; k < ; (21) 
and the combination of Eqs. (20) and (21) yields: 

-*=it£; < 22 » 

D-D.= (23) 

finally, the related distribution of metal abundance in long-lived stars can 
be determined following a standard procedure (e.g., Pagel & Patchett 1975; 
Caimmi 1981), as done in C01. In fact, all the considerations made in C01 
(Appendix C) are independent of the sign of the inhibition parameter, k, and 
Theorem 1 therein may be generalized in the following way. 

Theorem Given a simple model of chemical evolution with star formation 
inhibiting gas, long-lived (i.e. life time longer than the age of the system) 
lower stellar mass limit, m m f, and inhibition parameter, k (where positive 
and negative values correspond to star formation inhibiting and enhancing 
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gas, respectively) , it is equivalent to any model of the same kind and with 
same values of parameters, except the above mentioned two, (m m f) n and n n , 
which are defined by the relations: 

(m m f)i < (m mf ) n < m mf ; k > ; (24a) 
(mmf)n < (m m f)i < m mf ; -1 < k < ; (24b) 

1 + K n 
1 + K 

where (m m f)i is related to K n = 0, i.e 
enhancing gas. 

Aiming to an application of the above results to the chemical evolution of 
the Galactic spheroid, the value of the real normalized yield, p/O , the long- 
lived star and remnant mass fraction in a star generation, a, and the lower 
stellar mass limit, m m f, shall be taken from a fit to the EGD in the disk solar 
neighbourhood (COO) with regard to a solar oxygen abundance, O = 0.0056 

(coi). 

Two extreme values of the power-law IMF exponent, p, shall be consid- 
ered, namely: (i) p = 2.9, which is a fit to the IMF determined by Scalo 
(1986) or Miller & Scalo (1979), for m ~ m , and provides a good approx- 
imation also in terms of oxygen production (Wang & Silk 1993); in fact, a 
different model (with respect to the current one) for the chemical evolution 
of the Galactic halo also requires a Miller-Scalo IMF to fit the data (Lu et 
al. 2001), and (ii) p = 2.35, which coincides with the IMF determined by 
Salpeter (1955). Different fits to the EGD in the Galactic spheroid, with 
respect to the solar neighbourhood, would imply star formation inhibiting or 
enhancing gas. For deeper insight into the model, see COO, C01. 

Input parameters which remain fixed are listed in Tab. 7, where the 
indices, 2.9 and 2.35, denote the value of the power-law IMF exponent 
used in computing the corresponding quantities. The initial oxygen abun- 
dance assumed for field halo stars, (p Q H = 0.001, is consistent with a lower 
limit deduced from theoretical considerations, to allow pop. II star formation 
(Bromm & Loeb 2003). 

Different models can be obtained by changing a single remaining input 
parameter, which can be chosen as the normalized effective yield, p"/0 & , or 
the slope of the TGD, a, conform to Eqs. (14) and (15). 



= P - ; (25) 
p 

star formation neither inhibiting nor 
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Table 7: Values of input parameters of simple homogeneus models with star 
formation inhibiting or enhancing gas. The indices, 2.9 and 2.35, denote the 
value of the power-law IMF exponent used in computing the corresponding 
quantities. The indices, H and B, denote halo and bulge field star subsystem, 
respectively. 



p/Oq ■ 10 7 

{m m f) 2 .9 ■ 10 3 

(rnmfh.35 ■ 10 3 6 

"2.9 • 10 7 

"2.35 • 10 8 

Ho 1 

So 

D 

(PoH ■ 10 3 1 

4>oB o 

4>fH l 



3722 
4235 
9136 
3666 
9104 



20 




5 
6 



21 



Table 8: Values of parameters related to homogeneous simple models with 
star formation inhibiting or enhancing gas, corresponding to linear fits to the 
empirical differential metallicity distribution (EGD) in halo (H) and bulge 
(B) field stars, plotted in Fig. 4. The mean oxygen abundance (normalized 
to the solar value) of stars at the end of evolution is denoted as 0. Positive 
and negative k values mean star formation inhibiting and enhancing gas, 
respectively. Positive and negative D values mean star formation inhibiting 
gas and stars formed from additional gas with same composition as the pre- 
existing gas, respectively. Mass fractions are normalized to the initial mass. 





HI 




H2 




Bl 




B2 




p"/o 


1.3363 E- 


1 


9.2288 E- 


-2 


1.0857 E- 


-0 


8.1430 E- 


1 


a 


-3.2500 E- 


-0 


-4.7059 E- 


-0 


-4.0000 E- 


-1 


-5.3333 E- 


1 


b 


8.7771 E- 


1 


1.0396 E- 


-0 


1.1048 E- 


-1 


2.3465 E- 


■1 


K 


4.5169 E- 


-0 


6.9883 E- 


-0 


-3.2099 E- 


-1 


-9.4658 E- 


2 


H 


5.6657 E- 


4 


1.9899 E- 


-5 


7.5858 E- 


-3 


1.4905 E- 


3 


s f 


1.8116 E- 


•1 


1.2518 E- 


-1 


1.4616 E- 


-0 


1.1029 E- 





Df 


8.1828 E- 


•1 


8.7480 E- 


-1 


-4.6915 E- 


-1 


-1.0440 E- 


-1 





1.3406 E- 


1 


9.3268 E- 


-2 


1.2452 E- 


-0 


1.0064 E- 


■0 



A few cases are listed in Tab. 8, related to halo (H) and bulge (B) field 
stars, respectively, where represents the mean oxygen abundance (normal- 
ized to the solar value) of stars at the end of evolution. Positive and negative 
k values mean star formation inhibiting and enhancing gas, respectively. 
Positive and negative D values mean star formation inhibiting gas and stars 
formed from additional gas with same composition as the pre-existing gas, 
respectively. 

It is apparent that halo and bulge models demand star formation inhibit- 
ing and enhancing gas, respectively, to maintain (i) universal power-law IMF, 
and (ii) real normalized yield, p/0 & , unchanged with respect to a value de- 
duced from an acceptable fit to the EGD in the disk solar neighbourhood 
(COO; C01). More precisely, about 80-90% of the initial halo gas has to be 
inhibited from forming stars and, on the other hand, about 10-50% of the 
initial bulge mass has to be added as gas enhanced in forming stars, for 
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providing an acceptable fit to the related EGD. It results in a normalized ef- 
fective yield lower by about one order of magnitude in the halo, with respect 
to the bulge. 

The TGD deduced from models H1-H2 and B1-B2 is represented in Fig. 6 
and compared to the corresponding EGD in connection with halo (top panels) 
and bulge (bottom panels) field stars, both in presence (left panels) and in 
absence (right panels) of [O/Fe] plateau. 

With regard to the halo, homogeneous simple models in presence of star 
formation inhibiting gas, provide an acceptable fit with values of input pa- 
rameters as listed in Tab. 7 and 0.092 < p" /O < 0.134, both in presence 
and in absence of [O/Fe] plateau. 

With regard to the bulge, homogeneous simple models in presence of 
star formation enhancing gas, provide an acceptable fit with values of input 
parameters as listed in Tab. 7 and 0.81 < p"/0 & < 1.09, both in presence 
and in absence of [O/Fe] plateau. 

Further inspection of Fig. 6 shows that, in both cases, there is a slight 
deficiency (strenghtened in absence of [O/Fe] plateau), in the number of 
stars observed below a treshold, 0.004 or [Fe/H]f« —3.0 with regard to 
the halo, and « 0.5 or [Fe/H]« —0.6 with regard to the bulge. In other 
words, a G-dwarf problem seems to exist for both the halo (Hartwick 1976; 
Prantzos 2003) and the bulge (Ferreras et al. 2003). 

The TGD deduced from models H1-H2 and B1-B2 is represented in Fig 7 
and compared to the corresponding EGD in connection with the complete 
sample (N = 149) of globular clusters considered above (top panels), and a 
reduced sample (N = 107) with only old halo and bulge/disk objects retained 
(bottom panels), both in presence (left panels) and in absence (right panels) 
of [O/Fe] plateau. With regard to halo clusters, homogeneous simple models 
with star formation inhibiting gas, provide an acceptable fit leaving aside 
the occurrence of a G-dwarf problem. Accordingly, field halo stars and halo 
globular clusters underwent similar chemical evolution, which is consistent 
with recent results on the comparison of elemental abundance ratios (Pritzl, 
Venn & Irwin 2005). 

With regard to bulge/disk clusters, homogeneous simple models with star 
formation enhancing gas, provide an acceptable fit below a treshold in metal 
abundance, ps 0.6 or [Fe/H]» 0.5, while a G-dwarf problem appears when 
the treshold is exceeded. On the other hand, it is not the case for bulge field 
stars, which implies that the formation of bulge/disk metal-rich globular 
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clusters was inhibited (in the sense that proto-cluster stars turned into field 
stars) and/or tidal disruption took place. Accordingly, bulge (and possibly 
thick disk) field stars and bulge/disk globular clusters underwent similar 
chemical evolution. 

Clusters within a restricted metallicity range, -1.2 ~[Fe/H]~ -0.8, are 
fitted by both halo and bulge models, which allows to shift on the left the 
transition from halo to bulge/disk, from [Fe/H]f« —0.8 to [Fe/H] ~ —1-2, con- 
sistent with a restricted range in age shown by metal- rich clusters, [Fe/H]~ 
— 1.0, see Fig. 1. 

With regard to Galactic spheroid field stars, the differential version of 
Eq. (13a) reads: 



, , dN , dN H + dN B 
ib = log = log 

~N H dN H N B dN B ' 



log 



(26a) 



N N H d(f) N N B dt 
N H = N FH + (N gc )h ; N B = N FB + (N GC ) B ; (26b) 

which includes the contribution from globular clusters. The combination of 
Eqs. (10), (14), and (26) yields: 



ip = log 



M H M B 

— exp w (a H (f) + M + exp lo (a B + b B ) 



(27) 



where, in general, exp n x = x u and exp c x = e x according to the usual nota- 
tion. Then the TGD related to the Galactic spheroid is expressed by Eq. (27). 

Two alternatives shall be considered, Si and S 2 , according if the coeffi- 
cients, a and b, are taken from cases HI, Bl, and H2, B2, of Tab. 8, respec- 
tively, and M B /M = 0.1, M B /M = 0.9, are assumed to a good extent. The 
resulting TGD is plotted in Fig. 8 and compared to its empirical counterpart, 
represented in Fig. 5 (bottom panels), both in presence (left panels) and in 
absence (right panels) of [O/Fe] plateau. A discontinuity is exhibited by the 
TGD at = 0.20, where bulge formation is assumed to start. Though a G- 
dwarf problem seems to exist for both the halo and bulge, the current model 
provides a viable interpretation to the occurrence of three extremum points, 
two maxima and one minimum, in the EGD. 
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5 Inhomogeneous, simple models 



As a viable alternative to homogeneous simple models, let us take into con- 
sideration inhomogeneous simple models already used in previous approaches 
(COO; C01). In this view, the proto-halo/proto-bulge is conceived as struc- 
tured into a number of discrete, entirely gaseous, identical regions, and a 
subsystem of long-lived stars, remnants, and star formation inhibiting gas, 
which have been generated earlier. The evolution occurs via a sequence of 
identical time steps. At the beginning of a step, star formation stochasti- 
cally takes place in a subclass of "active" regions, as described by simple 
homogeneous models, while the others remain "quiescent". At the end of a 
step, high-mass stars have died whereas low-mass stars have survived upto 
today, according to instantaneous recycling approximation. In addition, the 
enriched gas left from active regions is instantaneously mixed with the un- 
enriched gas within quiescent regions, to form a new set of identical regions 
for the next step. For further details, see COO. 

The gas oxygen abundance, averaged over the whole system, still shows a 
monotonic increase with time, and depends on the gas mass fraction, accord- 
ing to Eq. (17), provided the fraction of active regions equals the probability 
of a region to be active (expected evolution; for further details, see COO). 
Taking into consideration also star formation inhibiting or enhancing gas, 
the effective yield takes the expression (C01): 



concerning an active region, p" is the effective yield due to star formation 
inhibiting or enhancing gas, expressed by Eq. (18), and p! R the gas mass 
fraction at the end of a step; concerning the system, q may be thought of as 
an effective gas mass fraction within a region at the end of a step, i.e. the 
mean gas mass fraction averaged on both active and quiescent regions, and 
X is the probability of a region being active at the beginning of a step. For 
further details, see COO. 

Though the TGD cannot be analytically expressed in the framework of 
inhomogeneous simple models, still it can be done in the special case of 
expected evolution, with regard to the starting point, tp(<f> ) = ipo, and the 





(28) 



g = l-x(l-A4) I 0<g<l ; 



(29) 
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ending point of the initial step, ip{<p + A0 O ) = ipi, where A0 O is the net 
oxygen abundance (normalized to the solar value) increase in gas component 
at the end of the first step. The result is (C01): 



= ^(0o) = log 



dN 



d</> 



log 



O 







(30) 



i>\ = ^(0o + A0 O ) = log 



1 — exp 




1 AN 
NA<p o 
A0 



V 







+ log 



XP-o 



fl - flf 



(31) 



where AN Q is the number of stars belonging to a selected spectral class, 
within the oxygen abundance range, O < <fi < <p + A0 O , and A0 O = A<f> 
independent of the step, in the case under discussion of expected evolution. 

An approximate expression of the TGD related to the first step, ipi, with 
the terms up to the second order retained, reads (C01): 



i>i = v(0 o +A0 o ) 



1 O A0 



In 10 p" 2 



flog 







O©A0 
p" 2 



«l ; 

(32) 

which is valid also in the general case, provided A0 is replaced by A0 O and 
the probability, x, by the relative frequency, u Q = k /n , being k Q and n the 
number of active and all regions, respectively, with regard to the first step. 

The values of some parameters related to the expected evolution, con- 
cerning cases HI, H2, Bl, B2, are listed in Tab. 9. The parameters appearing 
therein, which equal their counterparts corresponding to homogeneous mod- 
els (Tabs. 7 and 8), must be connected with active regions, and for this reason 
some corresponding output parameters show different values. For further de- 
tails, see Appendix A. New input parameters are: the effective yield, p', and 
the normalized oxygen abundance increase at the end of a step with regard 
to the whole system, A0, and to an active region, A<p' R , respectively. 

Concerning an active region, parameter values at the end of a step are 
calculated using a homogeneous simple model with star formation inhibiting 
or enhancing gas. Let fi' R , s' R , and D' R , be related to star formation allowing 
gas, long-lived stars, and star formation inhibiting or enhancing gas, from or 
enhanced in forming stars, respectively. Taking star formation inhibiting gas 
and leaving other parameters unchanged, makes no variation in the final gas 
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Table 9: Values of parameters related to the expected evolution of inhomo- 
geneous simple models, in connection with two different cases for the halo, 
HI and H2, and for the bulge, Bl and B2, respectively. The indices, 2.9 and 
2.35, denote values related to the power-law IMF exponent, p, in computing 
the corresponding quantities. With regard to the parameter, ipi, upper and 
lower values are calculated by use of Eqs. (31) and (32), respectively. The 
effective yield, p', is related to inhomogeneities in oxygen abundance due to 
the presence of active and quiescent regions, whereas oxygen is uniformly 
distributed within active regions. The lower part of the table is related to 
models with star formation inhibiting or enhancing gas. Parameters not re- 
ported therein have same value as in the upper part, with the exception of 
p, a, and m m f, which are listed in Tab. 7 together with other parameters not 
appearing here. The effective yield, p" , related to star formation inhibiting 
or enhancing gas, is listed as p in the upper part of the table. The effective 
yield, p', due to the presence of both star formation inhibiting or enhancing 
gas within active regions, and star formation precluding gas within quies- 
cent regions, is listed with the same notation in the upper part of the table. 
The index, R, denotes a generic active region. The mean oxygen abundance 
(normalized to the solar value) of stars at the end of evolution is denoted as 
0. Positive and negative n-D values correspond to star formation inhibiting 
and enhancing gas, respectively. Mass fractions are normalized to the initial 
mass. 
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mass fraction, /ij, as it is unrelevant if gas is frozen into long-lived stars or 
inhibited from forming stars. Accordingly, the final star plus star formation 
inhibiting or enhancing gas mass fraction, Sf + Df — 1 — fj,f, also remains 
unchanged. On the other hand, s' R + D' R = 1 — /j,' R within active regions at 
the end of a step, where s' R /D' R = Sf/Df necessarily holds. The combination 
of the above relations yields: 

s f = \^s' R ■ (33) 

Df = \^ D ' R ; (34) 

which are listed in the lower part of Tab. 9. 

The TGD related to cases HI, Bl, (full lines) and H2, B2, (dashed lines) 
is compared in Fig. 9 to the corresponding EGD with regard to halo (top 
panels) and bulge (bottom panels) field stars, both in presence (left panels) 
and in absence (right panels) of [O/Fe] plateau. 

With regard to the halo, inhomogeneous simple models with star forma- 
tion inhibiting gas, provide an acceptable fit using values of input parameters 
listed in Tab. 7 and 0.0033 < p'/0 Q < 0.0054, in presence and/or in absence 
of [O/Fe] plateau. 

With regard to the bulge, inhomogeneous simple models with star forma- 
tion enhancing gas, provide an acceptable fit using values of input parameters 
listed in Tab. 7 and 0.56 < p'/0 Q < 0.77, in presence and/or in absence of 
[O/Fe] plateau. 

In dealing with active regions, the normalized oxygen abundance at the 
end of evolution, <f>f, has to be replaced with the normalized oxygen abun- 
dance at the end of a step, fa + Afa R . 

Further inspection of Fig. 9 shows that, with respect to homogeneous 
simple models plotted in Fig. 6, the fit is more or less unchanged for the 
halo, and slightly improved for the bulge. In any case, a G-dwarf problem 
still remains. 

The TGD related to cases HI, Bl, (full lines) and H2, B2, (dashed lines) 
is compared in Fig. 10 to the corresponding EGD with regard to the whole 
sample (N = 149) of globular clusters considered above (top panels), and 
a reduced sample (N = 107) of old halo and bulge/disk objects (bottom 
panels), both in presence (left panels) and in absence (right panels) of [O/Fe] 



28 



plateau. The trend looks like its counterpart exhibited by homogeneous 
simple models, and similar considerations can be made. 

The TGD related to the Galactic spheroid has to be numerically com- 
puted using Eq. (26) together with HI, Bl; H2, B2; models listed in Tab. 9, 
yielding cases SI, S2, respectively. Fractional masses equal to M^jM = 0.1, 
Mb/M = 0.9, have been assumed to a good extent. A comparison with the 
related EGD (Fig. 5, bottom panels) is made in Fig. 11 both in presence (left 
panels) and in absence (right panels) of [O/Fe] plateau. The trend looks 
like its counterpart exhibited by homogeneous simple models, and similar 
considerations can be made. 

The temporal behaviour of (allowing star formation) gas oxygen abun- 
dance (TAMR), normalized to the solar value, log0 =[0/H], related to mod- 
els HI and Bl, are plotted in Figs. 12 and 13, in presence and in absence of 
[O/Fe] plateau, respectively. Also plotted therein are the data coming from a 
sample of 55 Galactic globular clusters (De Angeli et al. 2005) but expressed 
in terms of absolute ages (De Angeli 2005), with same captions as in Fig. 1. 
The conversion from [Fe/H] to [O/H] has been made using Eqs. (1) and (2) in 
presence and in absence of [O/Fe] plateau, respectively. It is worth recalling 
that simple models of chemical evolution used in the current paper are time 
independent, with regard to the TGD. Accordingly, the initial and the final 
time, together with the time step in the case of inhomogeneous models, can 
be selected for best fitting the empirical age-metallicity relation (EAMR). 

With regard to the halo, initial and final values, ([O/H], T/Gyr)= (—3, 
12.5) and (0, 8.0), respectively, have been chosen, together with a time step, 
AT/Gyr= 1.125, for a total of four. The related TAMR is represented by 
four full curves in Figs. 12 and 13. 

With regard to the bulge, initial and final values, ([O/H], T/Gyr) = 
(—0.70, 10.5) and (0.74, 10.0), respectively, have been chosen, together with 
a time step, AT/Gyr= 0.125, for a total of four. The related TAMR is rep- 
resented by a single full curve in Figs. 12 and 13, the remaining three being 
out of scale on the right. 

The comparison between EAMR (related to globular clusters) and TAMR 
(related to both globular clusters and field stars) shows that surviving (or 
formed) halo globular clusters belong to all model time steps, while surviving 
(or formed) bulge/disk globular clusters belong to the first step only and, in 
fact, are coeval within the age uncertainties. 

In the framework of inhomogeneous models, star formation allowing gas 
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has to be able of generating both active and quiescent regions. With regard 
to models H and B listed in Tab. 9, almost all star formation allowing gas has 
been exhausted at the end of evolution, leaving about 18-13% of long-lived 
stars and 82-87% of star formation inhibiting gas in the halo, and about 72- 
91% of long-lived stars from primeval gas and 28-9% from inflowing gas in 
the bulge, related to the final mass. The mean normalized oxygen abundance 
within stars is about 0.13-0.093 and 0.94-0.83 in the halo and the bulge, 
respectively. It is worth of note that the Galactic spheroid (halo and bulge) 
would be mostly luminous for a IMF exponent, p = 2.9, and mostly dark for 
p = 2.35, owing to different lower stellar mass limits, by more than an order 
of magnitude. On the other hand, the assumption of a universal IMF implies 
the same yield, and lower stellar mass limit, for the spheroid and the disk. 

The comparison between homogeneous (Tab. 8) and inhomogeneous 
(Tab. 9) simple models, shows no substantial change in TGD. It is worth 
emphasizing that the only difference is in the physical process related to the 
effective yield, which is due to star formation inhibiting or enhancing gas for 
homogeneous models and, in addition, to the presence of active and quiescent 
regions for inhomogeneous models. On the other hand, the related TAMR 
cannot provide an acceptable fit in the former alternative, but can do in the 
latter, as shown in Figs. 12 and 13. 

6 Discussion 

According to the results of Sections 3 and 4, halo and bulge EGD, taken sep- 
arately, are consistent with both homogeneous (Fig. 6) and inhomogeneous 
(Fig. 9) simple models of chemical evolution. The former alternative has been 
widely discussed (e.g., Hartwick 1976; Caimmi 1981, 1982; Ryan & Norris 
1991). The latter one allows a direct comparison with its counterpart related 
to the disk solar neighbourhood (COO), under the assumption of a universal, 
power-law IMF. 

As regards the disk solar neighbourhod, a lower stellar mass limit, ex- 
ceeding the theoretical Jeans stellar mass (0.007 < fhj < 0.01), occurs for 
a power-law exponent p = —2.9, which is a fit to the Scalo (1986) IMF for 
m ~ m , concerning both mass distribution and oxygen production (Wang 
& Silk 1993). A less steep Salpeter (1955) IMF, implying p = —2.35, is 
marginally consistent with the theoretical Jeans stellar mass (Tab. 7). On 
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Table 10: Comparison between parameters of inhomogeneous simple mod- 
els, related to halo and bulge (cases HI and Bl of Tab. 9) and disk solar 
neighbourhood (C01, Tab. 5). An universal power-law initial mass function 
(IMF) is assumed in all cases, which leaves other parameters i.e. p, a, and 
m m f, unchanged. Gas and star mass fractions, /if, Df, and Sf, are related 
to the initial mass with regard to the bulge, where total mass conservation 
is violated by gas inflow. 
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the other hand, the occurrence of stellar wind would reduce oxygen nucle- 
osynthesis by a factor of about 2.5 (e.g., Wang & Silk 1993) which, in turn, 
would raise the lower stellar mass limit (e.g., C01). 

The more relevant parameters of inhomogeneous simple models related 
to halo, bulge, and disk solar neighbourhood, are listed in Tab. 10. The 
related values are taken from Tab. 9 (models HI and Bl) and C01 (Tab. 5), 
respectively. An universal power-law IMF is assumed, <f>(rh) oc fnTp, which 
makes no change in value of the physical parameters, p, a, and m m f. 

With regard to active regions, the probability of star formation is x ~ 1 
for the halo, x < 1 f° r the bulge, and ^ < 1 for the disk. Accordingly, 
the gas mass fraction left at the end of a step, fj/ R , is close to zero for the 
halo, about one seventh for the bulge, and about one third for the disk. The 
effective yield, p", is about six times larger in the disk (where it coincides 
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with the real yield) and in the bulge (where star formation is enhanced by 
inflowing gas with same composition as in the pre-existing gas), than in the 
halo (where gas is partially inhibited from forming stars). 

With regard to the whole system, a similar trend occurs. The ratio of 
gas mass fraction at the end and at the beginning of a step is q <C 1 for the 
halo, q < 1 for the bulge, and q ~ 1 for the disk. Accordingly, the (allowing 
star formation) gas mass fraction left at the end of evolution, /if, is close to 
zero for the halo and about four thousandths for the bulge, and about one 
third for the disk. The effective yield, p', is about one half thousand times 
larger in the bulge (where it is due to both inhomogeneous star formation 
and gas inflow), and about one thousand times larger in the disk (where it 
is due to inhomogeneous star formation) than in the halo (where it is due 
to both inhomogeneous star formation and inhibition from forming stars). 
With regard to the total mass at the end of evolution, it is left about 18% of 
long-lived stars (including remnants) and 84% of gas inhibited from forming 
stars in the halo; about 72% of long-lived stars from primeval gas and 28% 
from inflowed gas in the bulge; about 71% of long-lived stars and 29% of gas 
allowing star formation in the disk. 

The TGD related to both homogeneous and inhomogeneous simple mod- 
els provides an acceptable fit to the EGD related to the Galactic spheroid 
(Figs. 8 and 11) and, in particular, a non monotonic trend is reproduced. 
While models assume that star formation at the beginning of halo and bulge 
evolution starts abruptly with constant efficiency, the data seem to indicate 
a somewhat gradual rate, with increasing efficiency In other words, the oc- 
currence of some physical process (not necessarily gas infall) seems to inhibit 
halo and bulge star formation in early times. Accordingly, the history of 
each Galactic subsystem (halo, bulge, disk) could be conceived as made of 
two distinct phases, namely (i) assembling, where star formation efficiency 
is gradually increasing, and (ii) stabilization, where star formation efficiency 
maintains constant. 

At the end of halo evolution, the fractional gas and star mass predicted 
by the model, are: 



fif + D f ; 



(35) 




= s f ; 



(36) 



32 



where Mho is the initial halo mass. The combination of Eqs. (35) and (36) 
yields: 

M Hgas = H + Df M HstaTS ; (37) 
s f 

where M Hstais is the current halo mass. 

At the end of bulge evolution, the fractional gas plus stars and star mass 
due to gas inflow, predicted by the model, are: 



M B stars + M B gas 



M 



Bo 



Hf + s f ; (38) 



^r SL = - D f ; ( 39 ) 

MBo 

where Mb is the initial bilge mass and M^ stars is the star mass due to gas 
inflow, at the end of evolution. The combination of Eqs. (38) and (39) yields: 

Ms staJS = ^^-(M Bstars + M Bgas ) ; (40) 

where M# st ars + MB gas is the current bulge mass. 

With regard to models Hl-2, Bl-2, the values of fractional masses, ///, 
Sf, Df, are listed in Tab. 9, and taking a current halo mass Mh = M# sta rs = 
10 9 m o and a current bulge mass Mh = Ms stars + Ms gas = lO lo m , it is 
found the following: 

M Hgas = 4.5-7.0 lO 9 m ; (41) 
^staxs = 3-2-0.95 lO 9 m ; (42) 

which is in agreement with the idea, that a fraction of the current bulge 
mass inflowed from the halo. The disk might have evolved separately for two 
orders of reasons. 

First, the empirical distribution of angular momentum in halo, bulge, 
thick disk, and thin disk stars seems to be consistent with a decoupled dy- 
namical evolution of the halo and the thick disk i.e. dissipative halo-bulge 
and thick disk-thin disk collapse (Wyse & Gilmore 1992; Ibata & Gilmore 
1995). Accordingly, the chemical evolution of the above mentioned subsys- 
tems could also have been decoupled. 
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Second, the disk mass (gas + stars) is estimated as Md = MD Stars + 
M Dgas « (5.0 + 0.8)10 10 m o = 5.8 lO lo m (e.g., Prantzos & Silk 1988), which 
cannot be provided by outflowing halo gas for more than about one tenth of 
the above value, according to Eqs. (41) and and (42). 

On the other hand, a continuous transition seems to exist from an ex- 
tended (R ~ 20kpc), pressure-supported halo, to an inner, flattened (R ~ 15 
kpc), rotation-supported halo (Chiba & Beers 2000). In addition, no corre- 
lation appears between mean rotational velocity and metal abundance for 
values below [Fe/H]^ —1.7 or 0.08, whereas exhibits a linear trend for 
larger values. The mean rotational velocity exceeds the statistical fluctua- 
tions related to metal-poor stars from [Fe/H] ~ —1.3 or <p « 0.2 on, and 
the relative abundance of thick disk stars is less than 1% for [Fe/H]« —1.7. 
For further details, see Chiba & Beers (2000). It is worth noticing that the 
range, —1.7 ~[Fe/H]~ —1.3, with regard to ZW metallicity scale, is related 
to a larger age spread in globular clusters, as shown in Fig. 1. 

The chemical abundance of thick disk stars suggests a similar history to 
those of metal-rich ([Fe/H]» —1.3) halo stars (Prochaska et al. 2000). In 
addition, the thick disk abundance patterns show excellent agreement with 
the chemical abundances observed in metal-poor bulge stars, suggesting the 
two populations were formed from the same gas reservoir at a common epoch. 
For further details, see Prochaska et al. (2000). 

To include the disk in a qualitative discussion, let us assume that a bary- 
onic dark halo exists, or an equivalent mass was lost during bulge and disk 
formation (C01). Accordingly, the following relations hold: 

Mi Mi M vis 



. I = H,B,D\ (43a) 

M M vis M ' ' ' y ! 

M vis = M H + M B + M D ; (43b) 

M = M vis + M uns ; (43c) 



and, in addition: 

M H + M uns _ M stars _ M B + M D _ M gas 



M M ' M M 



(44) 



where M v ; s is the visible Galactic mass, M uns the unseen (baryonic) mass, 
and Mh + M uus is the total halo mass, including the baryonic dark subsystem 
or the lost mass. 
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Using the above quoted values for Mh, Mb, Md, it is found Mh/M y - 1s = 
0.01449, M B /M vis = 0.14493, M D /M vis = 0.84058, and M vis = 6.9 lO lo m . 
On the other hand, the values of fractional masses, /if, Sf, and Df, pre- 
dicted by the model, must be related to the total mass, M. Accordingly, the 
following relation holds: 



^^ = 0.98551 ; (45) 

-'"vis 



via Eq. (44), and: 



M JT = »f + D f-> ( 46 ) 

as predicted by the model via Eq. (35), provided the initial halo mass coin- 
cides with the initial Galactic mass. If, in addition, bulge and (thick) disk 
formation took place from the same gas reservoir, as suggested by observa- 
tions (Prochaska et al. 2000), it may safely be assumed that the proto-halo 
was the common reservoir, which implies Eq. (44). 

The combination of Eqs. (43a), (44), and (46) yields: 

M B + M D M vis 

~M—1^ = ^ + Df; (47) 

where the mass fractions, /if and Df, are listed in Tab. 9 with regard to 
models HI and H2. Using Eqs. (43) and (45), the following mass ratios are 
determined: 

Mgl = 0.83078-0.88768 ; (48a) 
M 

= 0.16922-0.11232 : (48b) 

M 
Mm 

0.012040-0.012865 ; (48c) 



M 
Mb_ 

M 

Md _ 
M 
Mh_ 

M uns 



0.12040-0.12865 ; (48d) 
0.69834-0.74618 ; (48e) 
= 0.071150- 0.11454 ; (48f) 
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which shows that the unseen (baryonic) halo, or the total amount of gas lost, 
has to be as massive as about the bulge. 

In the light of the current model, it is apparent that about 17% -11% of the 
total mass gave no contribution to bulge and disk formation. If this material 
is related to stars with a mass below a treshold, m , then the mass ratio of 
long-lived stars above and below the treshold must necessarily equal the mass 
ratio, M H /M uns . The further restriction to a power-law IMF, <f>(rh) oc m~ p , 
implies the validity of the relation (C01): 

rh 2 o- p -m 2 -f M uns ' ™ m ' ^ 

where m mr is the mass of the oldest halo stars which are currently leaving the 
main sequence. For a model halo formation from 12.5 to 8.0 Gyr ago, m mr = 
0.9 m Q is expected to be consistent with the theory of stellar evolution. In 
fact, a value m mr = 0.87 m is deduced by linear interpolation from results 
related to metal-free stars (Marigo et al. 2001), and the above value has to 
be raised due to the observational lack of halo stars with zero metallicity. 

Taking rh mr = 0.9, rh m f from Tab. 7, and using Eq. (48f), the treshold 
mass, rh , may be deduced from Eq. (49). The result is: 

m = 0.25048-0.26178 ; p = 2.9 ; (50a) 
m = 0.089514-0.10027 ; p = 2.35 ; (50b) 

accordingly, stars with lower mass (for fixed power- law IMF exponent, p) 
either escaped detection up today, or became unbounded to the halo. 

If, on the other hand, the unseen baryonic halo is gaseous and unbound 
to the Galaxy, the power-law IMF exponent related to rh m f = m , leaving 
the remaining input parameters unchanged, reads: 

p = 2.8089-2.8206 ; (51a) 
p = 2.6025-2.6205 ; (51b) 

related to Eqs. (50a) and (50b), respectively. Accordingly, the values p = 2.9 
and p = 2.35 may be regarded as fiducial limiting values. 

An alternative explanation demands a different IMF in different Galactic 
subsystems. A minor change could only occur in the lower stellar mass limit, 
owing to a larger Jeans stellar mass (e.g., Larson 1998). Accordingly, a gas 
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amount prescribed to form stars in the range below the Jeans stellar mass, 
m m f < m < m , would follow a different fate, being lost from the system 
(e.g., Binney et al. 2001). The inhibition of star formation in the mass 
range under discussion, could be due, in addition, to less efficient cooling in 
metal-poor proto-stars. In fact, the pre-main sequence life time would be 
increased, to exceed the main sequence life time of massive stars. Finally, 
low-mass proto-stars would be destroyed by supernovae and the related gas, 
which is expected to be weakly bound to the proto-Galaxy, would be lost 
during bulge and disk formation. 

7 Conclusion 

The empirical differential oxygen abundance distribution (EGD) in the Galac- 
tic spheroid has been deduced from three different samples involving (i) 
268 K-giant bulge stars (Sadler et al. 1996), and (ii) 149 globular clusters 
(Mackey & van den Bergh 2005) for which the iron abundance distribution 
is known, in addition to previous results (Caimmi 2001) related to (iii) 372 
solar neighbourhood halo subdwarfs (Ryan & Norris 1991). To this aim, two 
alternative [0/H]-[Fe/H] dependences have been used, according to Eqs. (1) 
and (2), respectively. The data have been fitted, to an acceptable extent, by 
both homogeneous and inhomogeneous simple models of chemical evolution. 

Under the assumption of a universal initial mass function (IMF) and 
same value of the true yield as in the disk solar neighbourhood, inhibition of 
halo star formation (implying gas outflow) and enhanchement of bulge star 
formation (implying gas inflow) have been demanded for fitting the EGD. 
On the contrary, no such gas outflow or inflow was requested to reproduce 
the EGD in the disk solar neighbourhood (COO). A power-law IMF has been 
considered, (f>(rh) oc m~ p , within the range, 2.35 < p < 2.9, and special 
effort has been devoted to the limiting cases, p = 2.9, which is acceptably 
close to Scalo IMF for m ~ m , and (ii) an exponent p = 2.35, which is the 
Salpeter IMF. In any case, it has been inferred that a more refined model 
involving an initially increasing star formation efficiency (but not necessarily 
implying gas infall) while assembling Galactic subsystems, could provide a 
better agreement with the data. 

Homogeneous models have been recognized unable in fitting the empiri- 
cal age-metallicity relation (EAMR) with regard to a homogeneous sample 
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of globular clusters De Angeli et al. 2005; De Angeli 2005), which shows a 
non monotonic trend characterized by large dispersion. On the other hand, 
inhomogeneous models have been shown an acceptable fit, provided globu- 
lar cluster formation occurred through four different steps in the halo and 
through a single step in the bulge/disk, unless clusters of later generation 
were disrupted. 

With regard to gas outflow from the halo, acceptable models made the 
following predictions. If spheroid and disk component underwent distinct 
evolutions, then a non negligible fraction of the bulge mass (from about one 
third to about one tenth) outflowed from the halo, for assumed Mh = 10 9 m 
and Mb = 10 10 m . If, on the other hand, spheroid and disk component 
underwent a common evolution, then an unseen baryonic halo (or equivalent 
amount of gas lost by the Galaxy) has been shown to be needed, for assumed 
Mfl = 5.810 10 m Q . The mass of the unseen halo has been found to be of 
the same order as the bulge mass. In addition, the treshold star mass below 
which the halo is not detectable (or the stars are unbound to the Galaxy) has 
been calculated as m pa 0.25 m for IMF exponent p = 2.9, and m pa 0.1m Q 
for p = 2.35; conversely, p pa 2.8 for lower stellar mass limit m m f = 0.25 m , 
and p pa 2.6 for m m f = 0.1 m . 
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Appendix 



A Correspondence between homogeneous and 
inhomogeneous simple models 

With regard to a generic system, let /x G , ///; <f> , 0/; be initial and final values 
of fractional (allowing star formation) gas mass and oxygen abundance nor- 
malized to the solar value, respectively. Let chemical evolution be described 
using either homogeneous or inhomogeneous simple models. Finally, let the 
generalized yield, p", related to the system in the former alternative, coincide 
with its counterpart related to an active region in the latter alternative. 

With regard to a generic active region at the end of a step, Eq. (17) 
reduces to (C01): 



where A(fr' R = <fr' R f — (fr' Ro , jj,' Ro = 1, and ji' R j = p! R . The combination of 
Eqs. (3), (17), and (52), in the case under discussion yields: 



which allows the calculation of the fractional gas mass ratio, /i///U , predicted 
by homogeneous simple models, in terms of parameters related to inhomoge- 
neous simple models with equal values of initial and final normalized oxygen 
abundance, <f> , and (frf, respectively. 

With regard to homogeneous simple models, the mean normalized oxygen 
abundance in long-lived stars at the end of evolution is (C01, Appendix D): 



A0' 




ln/4 ; 



(52) 



(frf ~ <fro = ln(/i///x ) 
A(fr' R ln/z^ 



(53) 



<fr = 0(0/) = <fro + (0/ - 4>o) 



u(l — \nu) — 1 




(54) 



(1 — u) lnw 



and the combination of Eqs. (3), (17), and (54) yields: 



(frf = (fro + 




u\nu 



l-u 




(55) 



where u = fif in the special case, ji — 1. 
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With regard to inhomogeneous simple models, the mean normalized oxy- 
gen abundance in long-lived stars at the end of the £-th step is (C01, Ap- 
pendix D): 

</V = — e + —7z rr- -, — A0 fl ; (56) 

where z/j is the relative frequency of active regions at the i-th step. 

In the special case of expected evolution, vi = x, the first term on the 
right-hand side of Eq. (56) reduces to (C01, Appendix D): 



l^i=Q Villi 



p' qliaq 



1 



(£ + 1)^(1 -q) 



1 



7 m 



(57) 



where p' is the generalized yield related to the system, and q may be thought 
of as an effective gas mass fraction within a region at the end of a step i.e. 
the mean gas mass fraction averaged on both active and quiescent regions. 
The combination of Eqs. (28) and (57) yields: 



— <Po 



Ei= 



o [ 



(£ + l)q e (l-q) 



(58) 



in terms of the generalized yield, p", related to an active region. The substi- 
tution of Eqs. (52) and (58) into (56) produces: 



(fie = 0o + 



p 

G 







1 + 



/^ln/^ (£ + -q) 



1 e+i 



(59) 



at the end of the £-ih step. 

In the special case where the system reduces to a single region and the 
evolution to a single step, /i///i = /i' R , £ = 0, Eq. (59) reads: 



+ 



p" 



1 + 



/4 



(60) 



which coincides with Eq. (55), related to homogeneous simple models. 
The above results may be reduced to a single statement. 

Theorem Given (i) an inhomogeneous simple model of chemical evolution 
with assigned values of initial normalized oxygen abundance, (fi , generalized 



40 



yield related to an active region, p", and gas mass fraction within an active 
region at the end of a step, /i' R ; and (ii) a homogeneous simple model of 
chemical evolution with initial normalized oxygen abundance, (p Q , generalized 
yield related to the system, p" , and gas mass fraction at the end of evolution, 
/j,f — fi' R ; the mean oxygen abundance in long-lived stars at the end of the first 
step in the former alternative, equals its counterpart at the end of evolution 
in the latter alternative. 
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Figure 1: Empirical age-metallicity relation (EAMR) from a sample of 55 
globular clusters (De Angeli et al. 2005) in terms of absolute ages (De Angeli 
2005), with regard to ZW (top panels) and CG (middle panels) metallicity 
scales, and values from Harris catalogue (Mackey & van den Bergh 2005; 
bottom panels), with (left panels) and without (right panels) error bars. 
Morphological types: OH (old halo) - triangles; YH (young halo) - squares; 
BD (bulge/disk) - diamonds; NCG6366 - asterisk. 
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Figure 2: The empirical, differential metallicity distribution (EGD) in glob- 
ular clusters, plotted with regard to a complete sample (De Angeli et al. 
2005; iV = 55, top panels) and a reduced sample (N = 42, bottom panels) 
with YH clusters removed, both in presence (left panels, CG metallicity cali- 
bration) and in absence (right panels, ZW metallicity calibration), of [O/Fe] 
plateau, respectively. The dotted vertical line marks the transition from halo 
(OH, YH) to bulge/disk (BD) morphological type, [Fe/H]=-0.8. The dis- 
tribution appears to be bimodal, with the occurrence of two maxima, close 
to the beginning of evolution and to the transition from halo to bulge/disk 
morphological type, respectively. 
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Figure 3: The empirical, differential metallicity distribution (EGD) in glob- 
ular clusters, plotted with regard to a complete sample (Mackey & van den 
Bergh 2005; N = 149, top panels) and a reduced sample (N = 107, bottom 
panels) with only OH clusters retained, both in presence (left panels) and 
in absence (right panels) of [O/Fe] plateau, respectively. The dotted vertical 
line marks the transition from halo (OH, YH) to bulge/disk (BD) morpho- 
logical type, [Fe/H]= —0.8. The distribution appears to be bimodal, with 
the occurrence of two maxima, close to the beginning of evolution and to the 
transition from halo to bulge/disk morphological type, respectively. 
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Figure 4: The empirical, differential metallicity distribution (EGD) in solar 
neighbourhood field halo subdwarfs (top panels), plotted with regard to a 
sample of 372 objects (Ryan & Norris 1991), and in K-giant field bulge stars 
in Baade's Window (bottom panels), plotted with regard to a sample of 268 
objects (Sadler et al. 1996), both in presence (left panels) and in absence 
(right panels) of [O/Fe] plateau, respectively. The dotted vertical line marks 
the transition from halo (OH, YH) to bulge/disk (BD) morphological type, 
[Fe/H]=-0.8. 
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Figure 5: The empirical, differential metallicity distribution (EGD) in the 
light Galactic spheroid (globular clusters and field halo stars; top panels), and 
in the massive Galactic spheroid (globular clusters, field halo stars, and field 
bulge stars; bottom panels), both in presence (left panels) and in absence 
(right panels) of [O/Fe] plateau, respectively. The distribution has been 
deduced using Eqs. (11), (12), and (13), from a sample of 148 (top panels) 
or 149 (bottom panels) globular clusters (Mackey & van den Bergh 2005), a 
sample of 372 field halo subdwarfs (Ryan & Norris 1991), and a sample of 
268 K-giant field bulge stars (Sadler et al. 1996). The dotted vertical line 
marks the transition from halo (OH, YH) to bulge/disk (BD) morphological 
type, [Fe/H]=-0.8. 
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Figure 6: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in halo (top panels) and bulge (bottom 
panels) field stars, both in presence (left panels) and in absence (right panels) 
of [O/Fe] plateau, respectively. The straight lines correspond to homogeneous 
models, H1-H2 (top panels) and B1-B2 (bottom panels) defined in Tabs. 7 
and 8. Crosses represent the data and related uncertainties, as in Fig. 4. The 
dotted vertical line marks the transition from halo (OH, YH) to bulge/disk 
(BD) morphological type in globular clusters, [Fe/H]=— 0.8. 
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Figure 7: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in globular clusters of all morphologi- 
cal types (top panels) and with only old halo and bulge/disk objects retained 
(bottom panels), both in presence (left panels) and in absence (right panels) 
of [O/Fe] plateau, respectively. The straight lines correspond to homoge- 
neous models, H1-H2 (more inclined) and B1-B2 (less inclined), respectively, 
defined in Tabs. 7 and 8. Crosses represent the data and related uncertainties, 
as in Fig. 3. The dotted vertical line marks the transition from halo (OH, YH) 
to bulge/disk (BD) morphological type in globular clusters, [Fe/H]=— 0.8. 
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Figure 8: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in the Galactic spheroid (top panels) 
and zoomed for low oxygen abundance (bottom panels), both in presence 
(left panels) and in absence (right panels) of [O/Fe] plateau, respectively. 
The curves correspond to models, HI, Bl; H2, B2; defined in Tabs. 7 and 8 
and combined via Eq. (27). Crosses represent the combination of the data 
and related uncertainties via Eqs. (13), as in Fig. 5 (bottom panels). The 
dotted vertical line marks the transition from halo (OH, YH) to bulge/disk 
(BD) morphological type in globular clusters, [Fe/H]=— 0.8. 
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Figure 9: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in halo (top panels) and bulge (bottom 
panels) field stars, both in presence (left panels) and in absence (right panels) 
of [O/Fe] plateau, respectively. Full curves correspond to inhomogeneous 
models HI (top panels) and Bl (bottom panels), and dashed curves to H2 
(top panels) and B2 (bottom panels), defined in Tab. 9. Crosses represent the 
data and related uncertainties, as in Fig. 4. The dotted vertical line marks 
the transition from halo (OH, YH) to bulge/disk (BD) morphological type 
in globular clusters, [Fe/H]=— 0.8. 
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Figure 10: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in all sample globular clusters (top 
panels) and only old halo and bulge/disk objects (bottom panels), both in 
presence (left panels) and in absence (right panels) of [O/Fe] plateau, re- 
spectively. Full curves correspond to models HI and Bl, and dashed curves 
to H2 and B2, where H curves have larger slope with respect to B curves. 
Crosses represent the data and related uncertainties, as in Fig. 3. The dotted 
vertical line marks the transition from halo (OH, YH) to bulge/disk (BD) 
morphological type in globular clusters, [Fe/H]=— 0.8. 
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Figure 11: Comparison between theoretical (TGD) and empirical differential 
oxygen abundance distribution (EGD) in the Galactic spheroid (top panels) 
and zoomed for low normalized oxygen abundances (bottom panels), both 
in presence (left panels) and in absence (right panels) of [O/Fe] plateau, 
respectively. Full and dashed curves correspond to models HI, Bl; H2, B2; 
respectively, listed in Tabs. 7-9 and combined via Eq. (26). Crosses represent 
the combination of the data and related uncertainties via Eqs. (13), as in 
Fig. 5 (bottom panels). The dotted vertical line marks the transition from 
halo (OH, YH) to bulge/disk (BD) morphological type in globular clusters, 
[Fe/H]=-0.8. 
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Figure 12: Comparison between empirical (EAMR) and theoretical age- 
metallicity relation (TAMR) in presence of [O/Fe] plateau, according to 
Eq. (1). The data come from a sample of 55 globular clusters (De Angeli 
et al. 2005) but expressed in terms of absolute ages (De Angeli 2005). Other 
captions as in Fig. 1. Full and dashed curves are related to models HI and 
Bl, respectively. Halo star formation begins at ([O/H], T/Gyr) = (-3, 12.5) 
and ends at (0, 8.0), within four time steps, AT/Gyr=1.125. Bulge star 
formations begins at (-0.70, 10.5) and ends at (0.74, 10.0), within four time 
steps, AT/Gyr=0.125. The last three steps are out of scale on the right, 
and cannot be shown. The dotted vertical line marks the transition from 
halo (OH, YH) to bulge/disk (BD) morphological type in globular clusters, 
[Fe/H]=— 0.8. The dashed vertical line marks the minimum in the differen- 
tial oxygen abundance distribution (EGD) in the Galactic spheroid (Fig. 8), 
= 0.20 or log0 = [0/H]= -0.70. 




Figure 13: Comparison between empirical (EAMR) and theoretical age- 
metallicity relation (TAMR) in absence of [O/Fe] plateau, according to 
Eq. (2). The data come from a sample of 55 globular clusters (De Angeli 
et al. 2005) but expressed in terms of absolute ages (De Angeli 2005). Other 
captions as in Fig. 1. Full and dashed curves are related to models HI and 
Bl, respectively. Halo star formation begins at ([O/H], T/Gyr) = (—3, 12.5) 
and ends at (0, 8.0), within four time steps, AT/Gyr=1.125. Bulge star 
formations begins at (—0.70, 10.5) and ends at (0.74, 10.0), within four time 
steps, AT/Gyr=0.125. The last three steps are out of scale on the right, 
and cannot be shown. The dotted vertical line marks the transition from 
halo (OH, YH) to bulge/disk (BD) morphological type in globular clusters, 
[Fe/H]=— 0.8. The dashed vertical line marks the minimum in the differen- 
tial oxygen abundance distribution (EGD) in the Galactic spheroid (Fig. 8), 
= 0.20 or log0 =[0/H]= -0.70. 



